Spontaneous superconductivity and optical properties of high-T c cuprates 
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We argue that the high temperature superconductivity in cuprate compounds may be supported 
by interaction between copper-oxygen layers mediated by in-plane plasmons. The strength of the 
interaction is determined by the c-axis geometry and by the a6-plane optical properties. Without 
making reference to any particular in-plane mechanism of superconductivity, we show that the 
interlayer interaction favors spontaneous appearance of the superconductivity in the layers. At a 
qualitative level the model describes correctly the dependence of the transition temperature on the 
interlayer distance, and on the number of adjacent layers in multilayered homologous compounds. 
Moreover, the model has a potential to explain (i) a mismatch between the optimal doping levels 
for critical temperature and superconducting density and (ii) a universal scaling relation between 
the dc-conductivity, the superfluid density, and the superconducting transition temperature. 
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The layered structure of the cuprates is a well estab- 
lished faclsi related to a strong anisotropy in their elec- 
tronic and optical properties^. The common feature of 
all cuprates — famous for their ability to exhibit super- 
conductivity at high transition temperatures — is the 
presence of conducting CuC>2 layers separated by the so 
called charge reservoirs. These reservoirs are nearly in- 
sulating even in superconducting phased 

One of the most fascinating features of high tempera- 
ture superconductivity is a strong dependence of basic su- 
perconducting properties (including the most important 
quantity, the transition temperature T c ) on the c-axis 
structure of cuprates. In particular, there is a systematic 
dependence of the critical temperature on the number n 
of the closely packed Cu02 layers per a structural c-axis 
unit. In typical homologous series of superconducting 
cuprates the separation between the n multilayers lies 
in the range d ~ 6 — 15 A which is large compared to 
the spacing c; n t ~ 3.5 A separating individual layers in- 
side the multilayer. Since Ci n t is numerically close to the 
in-plane Cu-0 bond length, a(Cu-O) « 3.8 A, it seems 
reasonable to treat the multilayer as a single thick layer. 

Reflectance data indicate that the dielectric functions 
of the cuprates may qualitatively be described by the 
so-called two-fluid modelr*^ 
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which is a generalization of a simplest metallic dielectric 
function s(ui) = 1 — ujp/[u)(u> + *7p)]- The carriers are di- 
vided into normal and superfluid components which have 
different impacts onto optical and conducting properties 
of the cuprates. In Eq. uj ps and u> pn are the plasma 
frequencies of the superconducting and normal compo- 
nents, respectively, and £oo is the high-frequency limit 
of e. The relaxation processes of the normal-state elec- 
trons are described by (in general, frequency-dependent) 
scattering rate 7 pn (w). In ordinary metals the relax- 
ation rate is small compared to the plasma frequency 



(i.e., j^/ujp 1 « 5 x 10 -3 for aluminium). Contrary to 
the metals, both the finite conductivity and the relax- 
ation processes are essential for optical properties of the 
cuprates 2 " because in a typical cuprate 7pn(w P n)/w pn ~ 1. 

Both in-plane and out-plane reflectance data show a 
sharp drop at the frequencies higher than the plasma 
edge^ uip S I^Eoo, which determines a boundary of a trans- 
parency window. Despite existence of other descriptions 
of the optical conductivities^ we take Eq. as our 
starting formula for the sake of concreteness. 

The form of the dielectric function suggests the 
presence of plasmon-mcdiatcd phenomena at the energy 
scales governed by the characteristic plasma frequencies. 
These phenomena are usually studied with respect to 
the c-axis conductivity ("transverse plasmon")£& The 
importance of the a6-plasmons for a proper descrip- 
tion of the superconducting state in layered materi- 
als such as high-Tc cuprates was clearly stressed in 
Ref. IjJ In our complimentary study we show that de- 
spite the afr-plane plasmon is heavily damped^ it in- 
duces the spontaneous appearance of the superconduc- 
tivity in the layers. Philosophically, our approach re- 
sembles mechanisms based on the interlayer Josephson 
tunneling 8 and interplane Coulomb interaction^ as well 
as other approaches' including phenomenological mod- 
els of the Ginzburg-Landau type»ii 

The free energy per one d period per unit layer area 
S in the absence of external fields is given by a sum of 
the contributions from the normal (F n ) and the super- 
conducting (F s ) states of the layer, and the plasmon- 
mediated interaction between the multilayers (F p i): 

F = F n (uj pn ,j pn ) + F s (uj ps ) + F p i(uj ps ,u! pn ,j pn ) . (2) 

In each term we explicitly indicate the leading-order de- 
pendence on the optical parameters oj ps , to pn , and "f pn . 
Long-range modulations of the c-axis structure are ne- 
glected. We imply that the effect of the intralayer media 
is solely insulating 2 thus neglecting a small finite out- 
plane conductivity in the normal state. 
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The free energy of the normal state F n in © should 
depend on the optical parameters uj pn and 7 p „ related 
to a specific (in fact, model-dependent) behavior of elec- 
trons in individual CuC>2 layers. Since we would like to 
keep our approach as general as possible we exclude F n 
from our analysis concentrating on the difference in the 
free energies of the normal and superconducting states, 

SF = F a (u> ps ) + Fpi(wp S , uj pn , 7 pn ) - F p i(0, io m , j pn ) . (3) 

The free energy density of the superconducting state in 
Eq. J5J is written in the Ginzburg-Landau (GL) formic 
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where the integration is going over the volume Vi of the 
superconducting layer, 2m* is the effective mass of the su- 
perconducting carrier, and A, B are the GL phenomeno- 
logical parameters describing the behavior of the order 
parameter ip which is related to the density of the con- 
densed electrons n s — \ip\ 2 . The GL approach has known 
limitations, while being usually correct near the point of 
the superconducting transition. Universality arguments 
suggest that the GL parameters must depend on the in- 
trinsic layer properties while being generally less depen- 
dent on the intralayer structure. 

Under assumption of a spatial homogeneity of the or- 
der parameter ip and negligence of fluctuations of the 
electromagnetic field A, the supercurrent in Q vanishes 
and we arrive to the simple expression 
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where w n is the multilayer width. We used the relations 
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where Xl is the London penetration depth. 

Naively, if the layers were structureless very thick solid 
plates made of alike atoms interacting with the van der 
Waals potential U(r) = -kT 6 , then the interaction 
energy of the layers would be described by the well- 
known Hamaker formic, U p \(d) = —H/{12ird 2 ), where 
H = nir 2 p 2 is the Hamaker constant and p is the number 
density of atoms in the planes. 

None of the above assumptions is satisfied by the 
cuprate layers because of significance of retardation, re- 
laxation, dielectric absorbtion and geometrical suppres- 
sion effects. These effects are knowniS to diminish the in- 
teraction which still follows the Hamaker law U p \ cx dr 2 
Up to an inessential numerical coefficient 
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where e; n t is the intralayer dielectric function. The ge- 
ometrical factor G takes into account the "multilayer- 
insulator" periodic structure 2 ^ 
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FIG. 1: The suppression factors: geometrical, G, Eq. and 
dissipative, D, Eq. (lll)t . and their asymptotics at r — > 0, oo. 



where ip^ is the first derivative of the digamma function. 
The geometric factor G (shown in Fig.GJby the solid line) 
is a monotonically increasing function of r = w n /d. 

The interaction energy J7J is of an electromagnetic ori- 
gin. In an idealized limit of perfectly conducting layers 
the interaction energy may be imagined as the Casimir 
energy^ of the electromagnetic field stored between the 
layersi£. In the case of real materials we follow Refill 
and interpret Eq. @ as the interaction energy between 
the layers caused by the inlayer plasmons. The plasmons 
give a dominant contribution to the interlayer energy at 
short interlayer separations^ uj pn d <C 1. This condition 
is satisfied by typical cuprates (e.g., u pn d < 10~ 3 for the 
La2- a; Sr a ;Cu04 compound discussed below). The charac- 
teristic frequency^ £7 in Eq. (JJJ gives account of the ab- 
sorbtion spectra of the layers which, in turn, characterize 
the strength of the interaction between the afe-plasmons. 

The frequency f2 can be expressed via the dielectric 
function s{uS) evaluated at the imaginary axis oj = i£. 
The dispersion relation 12 expresses e(z£) via the conduc- 
tivity Rccr = lu Im. e / (4tt) at the real axis: 
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Thus, the interlayer interaction J7J is fixed by the dissipa- 
tive part of the in-plane conductivity, Re a, the interlayer 
dielectric parameter £i nt , and the c-axis geometry. 

Despite the scattering rate 7 pn in cuprates is a com- 
plicated frequency-dependent function 2 ^ one can approxi- 
mately evaluate the order of the characteristic frequency 
n in the normal state (u ps = 0) assuming that j pn is 
w-independent: Q = uj pn D(-y pn /tu pn ). The dissipative 
d£, (7) suppression factor [derived from Eqs. (JTJ, J7J), JSJl], 
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is a monotonically decreasing function of r = 7 pn /o 
We plot D(r) in Fig.QJby the dashed line. 
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To estimate the energy scales related to the plas- 
mon-mediated interactions, we consider L^-zSr^CuC^ 
(La214) compound. In La214 the CuC>2 layers are per- 
fectly flat and are separated by two LaO layers at the 
distance d — cq/2 6.6 A. The basic cell has the 
tetragonal structure clq x bo x cq with the base parame- 
ters ao « bo w 3.8 A. The condensation energy of the 
optimally doped (a; = 0.16, T c = 38 K) La214 com- 
pound is known to be 8 w 13 fieV per one atom 
of copper^ The normal state of the slightly under- 
doped La214 is characterized^ by the plasma frequency 
w pn ~ 6000 cm -1 while the typical scattering rate is of 
the order 7 P „ ~ 2000 cm -1 for frequencies higher than 
the a6-plane "pseudogap" u a b ~ 700cm~ 1 i 1 7i 1 ? Fixing 
the scattering rate to be constant and taking into ac- 
count the dissipative suppression factor, Z)(l/3) w 0.4, 
we get fiLa2i4 ~ 2500 cm" 1 . The characteristic fre- 
quency is of the order of a typical superconducting gap 
A ~ 10 — 50meV in the cupratesS, hQ^^u ~ 50meV. 

In order to avoid suspicious fine-tuning of parameters 
in Eq. (J7J) we roughly set e; n t ~ 1, cq ~ d (then ao ~ d 
and the geometrical suppression is G(l) sa 0.6). Then the 
plasmon-mediated interaction energy per copper atom is 
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This value is by an order of magnitude higher than the 



condensation energy U^f u ^ ~ lO^cond- °t ncr words, 
the condensation energy may well be explained by a 10% 
deviation in the interplane interaction, which in turn, 
should be related to a change of similar magnitude in 
the optical parameters of a cuprate as it cools down from 
critical to lower temperatures. In fact, the optical char- 
acteristics of cuprate compounds vary essentially in this 
ranged exhibiting, e.g., a sharp drop of the scattering rate 
7 ps and dominance of the superconducting component 
u>p s at T < T c . This argument stresses importance of the 
afe-plasmon mediated inter-layer interactions. Below we 
ignore all interlayer interactions except for Eq. J7J. 

The crucial property of the interlayer interaction term 
Upi is that it favors appearance of the scatterless su- 
perconducting component with uj ps ^ 0. To illustrate 
this property we expand the characteristic frequency O, 
Eq. |0I, at T — T c in powers of uj ps /uj pn <C 1, 

Q(u; ps ) = Lo p Ju + u 2 ■ (^Y + o\(^) 4 } }, (12) 
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The dimcnsionlcss coefficients certain functionals 

of the scattering rate % n (y) = l P n{y • u pn )/u> pn , e.g. 
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As one can see from Eq. <|13|) the second coefficient of 
the expansion (|12fl is always positive, u 2 > 0, regard- 
less of particularities of the scattering rate "f pn (uj). The 
behavior of at j pn = const is illustrated in Fig. [21 




FIG. 2: The characteristic frequency J7J fl vs the supercon- 
ducting frequency uj vs for various scattering rates j pn . 



Since the characteristic frequency f2 enters Eq. Q with 
the minus sign, the interaction energy J7J may provoke 
a tachyonic instability against emergence of the super- 
conducting condensate (JBJ, \yj\ 2 cx ui ps . In other words, 
the interlayer interaction JJJ supports the appearance of 
a superconducting (ip ^ 0) state provided the layers are 
intrinsically able — via any microscopic mechanism — 
to generate this superconducting state. 

At high temperatures the energy density associated 
with the superconducting condensate J2J, JSJ) is higher 
compared to the gain in free energy © which would be 
achieved by the plasmon interaction J7J. This makes the 
superconductivity energetically unfavorable. As the tem- 
perature decreases the coefficient A(T) in the GL free 
energy (J5J gets gradually smaller and at the certain tem- 
perature the quadratic u> ps term of the GL free energy (JSJ 
cancels the same term in the plasmon-mediated interac- 
tion J7J, (|12|) . This cancelation marks the critical tem- 
perature T c . At lower temperatures, T < T c , the overall 
coefficient in front of the quadratic term turns negative 
and the system becomes unstable against spontaneous 
development of the superconductivity u> ps ~ \tp\ ^ 0. 

The relation between the critical temperature (T c ), op- 
tical (jpn, w pn ) and geometrical (w n , d) parameters of the 
cuprates can be derived from Eqs. |0, J7J, (|12|) : 
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where the left hand side (LHS) contains optical and mi- 
croscopic parameters at T = T c while the right hand side 
(RHS) is of the purely geometrical origin. Below we list 
a few universal features of the cuprate superconductors 
which are described by Eq. (|14fl . 

Transition temperature T c vs n. The GL coefficient 
A(T) is a monotonically increasing function of temper- 
ature. Thus, the higher (lower) value of the RHS in 
(UU), the higher (lower) value of T c A The RHS of 
Eq. (|14fl is a linearly-rising function of w n at w n <C d. At 
w n = 0.4448 d the RHS has a maximum and then it de- 
creases as 1 /w n at w n ^> d (we used Eq. |(HJ| as well as the 
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asymptotics of G given in Fig.^). In n-layered cuprates 
the width of the multilayer is a monotonically rising 
function of n, which can approximately be estimated as 
w n = (n — l)(ci n t + 5 g ) + S c , where S g is the geometrical 
width of a single layer, ranging from <5 ff (La-214) = OA 
and <5 9 (Bi-2212) = 0.013A to <5 a (YBCO) = 0.274A, 19 
Cint ~ 3.5 A is the interlayer spacing inside the multilayer, 
and 8 C is a "coherence width" of external layers which 
should be of the order of the c-axis coherence length £ c 
(a few A). For typical crystallographic parameters the 
RHS of 1|14|) . and, consequently, the transition tempera- 
ture T c , are peaked around n max = 3. This behavior is 
in fact a universal feature of the homologous seriesi 

Transition temperature T c vs d. The r.h.s of (|14|) is 
a monotonically decreasing function of the separation d 
between the multilayers provided the other geometrical 
parameters of the c-axis structure are fixed. Thus, the 
larger d the lower temperature must be. This is another 
universal behavior observed in the cuprates^ 

Transition temperature T c vs x. One may expect 
that the highest T c is achieved at the doping x at 
which the density n s of the superconducting carriers 
is highest. However, this expectation is not confirmed 
experimentally^! the optimal doping for the transition 
temperature is noticeably lower compared to the one 
for the carriers (i.e., in La-214, Y-123, Bi-2212 cuprates 
one has x^° t w 0.16 < x"° t w 0.19). The plasmon- 
mediated interaction may explain this behavior. If the 
RHS of Eq. 114fl were independent of x then the maxi- 
mum temperature would be achieved at a certain value 
of T = r c (x"p t ) corresponding to the highest carrier den- 
sity. However, the interlayer distance d increases with 
the doping x^ lowering the plasmon interaction energy 
(proportional to the RHS of l|14|l). Thus, the equal- 
ity (|14|l is achieved at a lower value of the GL parameter, 
A[T c (x^ t )} < A[T c (x% t )], implying T c (x n o ; t ) < T c (x% t ). 

Scaling between T C) uj ps and dc- conductivity. At suf- 
ficiently low temperatures the normal component is al- 
most invisible in the dielectric functioniSiii In this case 
ft = irups/ (4\/2) [we used (fTTTft as well as the ^y pn — > 
asymptotic, Fig. ^ , and Q becomes linear in 



ducting frequency at T = as a minimum of ® : 
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Neglecting the quartic term in JSJ we get the supercon- 
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where we used Eq. I|14|l and disregarded the variation of 
the crystallographic parameters in the range of tempera- 
tures between T c and T = 0. The relation (|16|l links the 
ratio of the GL layer's parameter A &tT = T c and T = 
with both superconducting and normal optical proper- 
ties of the afr-planes. Note that (i) the LHS (RHS) of 
Eq. flfif depends solely on T — (T = T c ) quantities; 
(ii) the relation Hlofl is universal: it is does not depend 
on the c-axis structure and should hold for all cuprate 
materials with the same in-plane parameters. 

Since the relation (|l(j[) is dependent on Lo pn (T c ,uj) fur- 
ther analytical calculations are difficult. We notice, how- 
ever, that the integral (|13|l is saturated at low frequencies 
relevant to the normal-state dc-conductivity: 
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Therefore we take "f pn to be equal to its low-frequency 
extrapolation, and work in the "dirty" limit 7 P „ 3> Lu pn , 
arriving to U2 = Jpn/^pn- The requirement of the 
dirty limit is rather mild since even at j pn — uj pn the 
above relation holds within 10%. Next, we take the 
standard GL-like prescription A(T) = a(T + T c ), with 
< T c -C T c , where a and T c are the GL parameters de- 
scribing the free energy associated with the condensation. 
We consider the most energetically unfavorable case: the 
positive T c indicates that the layers alone are not able 
to support the superconductivity. Curiously, if the in- 
trinsic layer properties are related, T c = /3w ps (0) with 
(3 w 10 K • cm, then from (|l(jfl and (|17f) we get the scaling 
relation w ps (0) ~ 120<7d c 2^ observed experimentally^ 
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